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We investigate in detail the asymptotic properties of tachyon cosmology for a generic class of self- 
interaction potentials. The present approach relies on an existing (formal) mathematical equivalence 
between a Dirac-Born-Infeld (DBI) model - with a particular choice of the warp factor and of the 
potential for the DBI field - and standard tachyon cosmology, under an appropriate transformation 
of the DBI field. The above mathematical equivalence is used to generalize the dynamical systems 
study of tachyon cosmology to a wider class of self-interaction potentials beyond the (inverse) square- 
law and power-law ones. It is revealed that independent of the functional form of the potential, the 
matter-dominated solution and the ultra-relativistic (also matter-dominated) solution, are always 
associated with equilibrium points in the phase space of the tachyon models. The latter is always the 
past attractor, while the former is a saddle critical point. For inverse power-law potentials V oc <j)~ 2X 
the late-time attractor is always the de Sitter solution, while for sinh-like potentials V oc sinh _Q (\(j>), 
depending on the region of parameter space, the late-time attractor can be either the inflationary 
tachyon-dominated solution or the matter-scaling (also inflationary) phase. In general, for most 
part of known quintessential potentials, the late-time dynamics will be associated either with de 
Sitter inflation, or with matter-scaling, or scalar field-dominated solutions. 
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I. INTRODUCTION 



£ 



(i) 



Inflationary models of the universe have been studied 
from the string theory perspective because inflation pro- 
vides an explanation for the homogeneity and isotropy 
of the early universe. Additionally, recent astrophysical 
observations indicate us that the universe is presently un- 
dergoing a phase of accelerated expansion that has been 
attributed to a peculiar kind of source of the Einstein's 
field equations acknowledged as "dark energy" [lj. 1 The 
crucial feature of the dark energy which ensures an ac- 
celerated expansion of the universe is that it breaks the 
strong energy condition. The tachyon field arising in the 
context of string theory [4| provides one example of mat- 
ter which does the job. The tachyon has been intensively 
studied during the last few years also in application to 
cosmology [Ej]-14|; in this case one usually takes Sen's 
effective Lagrangian density for granted [|J : 
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and studies its cosmological consequences without worry- 
ing about the string-theoretical origin of the action itself. 
In the above equation <j> is the scalar tachyon field, VMts 
self-interaction potential, and (V0) 2 = d^4>d u 4>. 

Recently generalized nonlinear scalar field theories of 
the Dirac-Born-Infeld (DBI) type have been proposed 
[lj|-[l7j]- As it is the case for Sen's tachyon field, these 
theories have attracted much attention in recent years 
due to their role in inflation based on string theory [l8| . 
In the above scenarios the inflaton is identified with the 
position of a mobile D3-brane, moving on a compact 6- 
dimensional submanifold of spacetime (for reviews and 
references see [19|), which means that the inflaton is in- 
terpreted as an open string mode. The effects of such 
DBI-motivated fields in a Friedmann-Robertson- Walker 
(FRW) cosmology have been already studied by means 
of the dynamical systems tools, yielding scaling solutions 
when the equation of state of the perfect fluid is negative 
and in the ultra-relativistic limit Non-linear Born- 
Infeld scalar fields with negative potential have been also 
investigated within FRW cosmology 20]. The cosmo- 
logical dynamics of a DBI field when the potential and 
the brane tension are arbitrary power-law or exponential 
functions of the DBI field, has been studied very recently 
in [H|. 

A dynamical systems study of the FRW cosmology 
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within the simpler theories based on the tachyon La- 
grangian (H|) can be found in Ref.[l4j]. However the au- 
thors of [lj] were able to study self-interaction potentials 
of the power-law type only. For more general potentials 
the corresponding system of ordinary differential equa- 
tions in the phase space is not a closed system of equa- 
tions any more, and one has to rely on the notion of 
"instantaneous critical points" whose physical relevance 
is unclear. 

In the present paper we aim at studying the cosmolog- 
ical dynamics of the tachyon model given by ([TJ), for a 
wider variety of tachyon self-interactions potentials than 
in Ref . [3] • To this end we shall exploit a formal mathe- 
matical equivalence existing between a DBI field - hold- 
ing a particular relationship between the DBI potential 
and the warp factor - and the tachyon field, in combina- 
tion with the application of an approach formerly used, 
for instance in Ref.(22[, that allows to study a vast vari- 
ety of self-interaction potentials. The above mathemat- 
ical equivalence enables to extend the results of dynam- 
ical systems studies of DBI models - including results 
of studies already existing in the bibliography [14j - to 
tachyon field cosmology. We will be able, in particular, to 
note an equivalence between the results of Ref. [T3| within 
the context of tachyon cosmology, and those in [2l| for a 
DBI field with exponential potential and brane tension, 
not reported previously. 

The paper has been organized in the following manner. 
The mathematical aspects of the particular DBI cosmo- 
logical model considered, are given in section II. Section 
III is devoted to the study of the asymptotic properties 
of the above model for a generic variety of self-interaction 
potentials, through the application of the dynamical sys- 
tems tools. In section IV we discuss the relevant aspects 
of the mathematical equivalence among the DBI model 
whose dynamics has been discussed in the previous sec- 
tion, and the corresponding tachyon model, so that the 
results of section III can be safely translated to the case of 
standard tachyon cosmology. It is demonstrated, in par- 
ticular, that former studies within the standard tachyon 
model, constrained to power- law potentials only [3], can 
be generalized to exponential type of potentials. A de- 
tailed discussion of the main results of the paper is pre- 
sented in section V, while the conclusions are given in 
section VI. For self-consistency and completeness, an ap- 
pendix with the basic recipes of the application of the 
dynamical systems tools to cosmology has been added. 
Here we use natural units c = 8nG = 1. 



sdbi = - J d^^Mir 1 wvi + /W(v#) 2 

-f-\d>) + V{ct>)}, (2) 

where /(0) is the inverse of the brane tension (also ac- 
knowledged as the warp factor of the warped throat ge- 
ometry in the internal space) and, V(<p) is the poten- 
tial for the DBI field, arising from quantum interactions 
between the D3-brane associated with cf>, and other D- 
branes. In the bibliography one usually encounters given 
forms of the warp factor /(</>). For instance, f(<fi) = A0~ 4 
(A-constant), or, also f{4>) — const. As we will immedi- 
ately see, in the present paper we choose f(4>) to be re- 
lated with the potential of the DBI field in the following 
particular form: /(</>)• V(tf>) = 1. The above choice leads 
to significant simplification of the field equations. 

For a spatially fiat FRW metric and arbitrary f(c/>) and 
V(4>), the equation of motion of the DBI scalar field, cj>, 
can be written in the following way: 
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where H = a/a is the Hubble expansion parameter, the 
dot accounts for derivative in respect to the cosmic time, 
and the "Lorentz boost" 7l is defined in the following 
way: 
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Alternatively the equation of motion of the DBI field 
can be written in the form of a continuity equation: 



p (i> + 3H(p<p +P4,) = 0, 



(5) 



where we have defined the following energy density and 
pressure of the DBI scalar field: 



f 



7l - 1 



V(<f>). (6) 



II. BASIC EQUATIONS AND SET-UP 



Significant simplification of the above equations can be 
achieved after considering the following particular choice 
of the warp factor /(</>), 

f{4>) = i/v(4>), 

leading to a simplified form of the effective action J2|): 



Consider the following effective action for a DBI field, 



St = 



d*xy/\g\V(<l>)y/l + (y<t>)*/V(4>). (7) 
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The Einstein's field equations in (flat) FRW metric, 
sourced by a mixture of a perfect barotropic fluid with 
energy density and pressure p m and p m — co m Pm respec- 
tively (u! m is the equation of state parameter of the per- 
fect fluid which, for dust, vanishes), and of a non- linear 
DBI scalar field, that can be derived from the action (0 , 
are the following: 



3H — p m + p^, 

2H = - (u> m + 1) Pm - (P4> + P<j>), 

Pm = -3 (LO m + 1) Hp m , 

4> + 3j- 2 H<p = -d^V{l - Ztf/W), 

where 

P4> = J v > P<t> = -7 -1 V, 
and the modified Lorentz boost 7 is defined as: 



(8) 



1 - ft/v 



(9) 



From now on we shall call the model described by equa- 
tions (|g|9"|i as "modified tachyon cosmology" (MTC). 

It is remarkable that the Lagrangian densities (0 and 
that in the action (0 are equivalent under the change of 
field variables 



do 



vnf) 



(10) 



Actually, under the transformation ()10|) the equations 
(0 above transform into the cosmological equations of 
standard tachyon cosmology (STC), in particular: 



2H = - (bJ m + 1) p rn - 



l~tp 2 



d v V 
V 



(11) 



As we shall show in the following sections, the above 
change of variable will allow us to translate the results 
of the dynamical systems study of a DBI field with the 
particular choice f(<fr) • V(4>) = 1 (next section), to stan- 
dard tachyon cosmology given by the Lagrangian density 
(0, thus allowing for a generalization of previous studies 
- see [lij for a dynamical systems study of power-law 
tachyon potentials - to a wider class of self-interaction 
potentials beyond the inverse square-law and power-law 
ones. 



III. DYNAMICS OF A DIRAC-BORN-INFELD 
FIELD (CASE f{<j>) • V((/>) = 1) 

Finding exact solutions of the cosmological equations 
|S]) is, in general, a very difficult task. That is why we 



will rely on the dynamical systems tools to investigate 
the asymptotic structure of the DBI/tachyon cosmology. 
To this end we will apply the concise recipes given in 
the appendix (section IVII|) . The goal will be to write 
the system of cosmological equations (JSJ in the form of 
an autonomous system of ordinary (ODE) as described 
in the appendix, so that one could associate such im- 
portant dynamical systems concepts as past and future 
attractors (also saddle equilibrium points), with dynam- 
ical configurations - solutions - of the models. This is 
a powerful approach to uncover the most generic classes 
of solutions that are allowed by them. In order to build 
an autonomous system out of the system of cosmological 
equations (|8l9p we introduce the following dimensionless 
phase space variables: 



V 



, y 



7 = 



(12) 



After this choice of phase space variables we can write 
the following autonomous system of ordinary differential 
equations: 



(x 2 

y 
2 



l){3x + V3(d^\nV)y}, 
W%xy{d < j>\sxV) + 37 m + 
Zy 2 (x 2 - j m ) 



(13) 



(14) 



where we have introduced the barotropic index for mat- 
ter 7 m = LJ m + 1 (do not confound with either 7^ or 7 
which play the role of the standard and modified Lorentz 
boosts respectively), and the tilde denotes derivative with 
respect to the time variable r = In a — properly the num- 
ber of e-foldings of expansion. While deriving the ordi- 
nary differential equations ([T5|). and ([H]). the following 
Friedmann constraint has also been considered: 



3iF 



(15) 



It will be helpful to have other parameters of observa- 
tional importance such as fl^ = p$ /3H 2 - the scalar field 
dimensionless energy density parameter, and the equa- 
tion of state (EOS) parameter lu^ — -p^j p^^ written in 
terms of the variables of phase space: 



VT 



1. 



Additionally, the deceleration parameter q 
H'/H): 



7m + 



y 2 {x 2 - 7m) 

VT^ 2 



(16) 
-(1 + 

(17) 
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A. Exponential Potential 

For an exponential self-interaction potential of the 
form: 

since d<p In V = —A = const, then the equations (|13p and 
(P3J form a closed autonomous system of ODE: 

x' = (x 2 - l){3x- \/3Ay}, 

?7 = g | 3 Tm - V3Axy + — l_ x2 J ■ ( 18 ) 

The phase space where to look for equilibrium points 
of the system of ODE (|T8]) - corresponding to the MTC 
model described by © - can be defined as follows: 

* = {(x,y) : -1 < x < 1, < y 4 < 1- x 2 }. (19) 

The equations (JTSJ) coincide with Eqs. (8,9) of [IJ], 
for a standard tachyon field with inverse square-law self- 
interaction potential of the form: 

V((j>) = V o 0- 2 . 

Hence, the same critical points as in [3] are found in 
the present case, this time for the exponential poten- 
tial. 2 This coincidence is easily explained through the 
equivalence between the present MTC model and stan- 
dard tachyon cosmology (STC) model, under the change 
of variable (TTU1) . Actually, as it will be shown in section 
IV, under (JTOj): 

V{4>) = V e~^ -> V(ip) =V (<p- ^T 2 . 

In the FigureQ]we show the trajectories in phase space 
for different sets of initial conditions for the model driven 
by an exponential potential. The free parameters have 
been arbitrarily set to 7 m = 0.25 and A = 1. Due 
to this choice of the parameters, the equilibrium points 
U = (±1,0) represent inflationary critical points - past 
attractors - in the phase space, while the matter-scaling 
solution is the late-time attractor. Points U are associ- 
ated with ultra-relativistic behavior since x — > ±1 => 
7 — > oo. 

Due to our definition of the variable x = <j)/Vv, the 
past attractor represents a scaling of the potential and of 
the kinetic energy of the tachyon scalar 

x = ±1 =► 4> 2 = v{4>). 



2 Recall, however, that in the above reference the variable x is 
defined in a simpler way: x = cf>. 




0- 
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FIG. 1: Trajectories in phase space for given initial data, for 
the exponential potential. The free parameters have been 
arbitrarily set to 7 m = 0.25 and A = 1. The matter-scaling 
solution is the late-time attractor, while the potential/kinetic 
energy scaling solution is the past attractor in the phase space. 



This means that, in case the inflationary past attractor 
be identified with early-time inflation, it can not be as- 
sociated with the slow-roll approximation which implies 
that the potential energy of the scalar field dominates 
over the kinetic one (the latter may be disregarded). 

In the next section we will focus on the asymptotic 
properties of MTC models driven by several others self- 
interaction potentials of cosmological relevance. 

B. Self-interaction Potentials beyond the 
Exponential one 

As long as one considers just constant and exponen- 
tial self-interaction potentials (8$ In V = and In V — 
const respectively), the equations (| 13[) and (TT41 form a 
closed autonomous system of ODE. However, if one wants 
to go further to consider a wider class of self-interaction 
potentials beyond the exponential one, the system of 
ODE (|13I14|) is not a closed system of equations any 
more, since, in general, d^, In V is a function of the scalar 
field cf>. A way out of this difficulty can be based on the 
method developed in [22j]. In order to be able to study 
arbitrary self-interaction potentials one needs to consider 
one more variable v, that is related with the derivative 
of the self-interaction potential through the following ex- 
pression 

v = -d^V/V = -dj In V. (20) 
Hence, an extra equation 

v' = -V3xyv 2 (T- 1), (21) 

has to be added to the above autonomous system of equa- 
tions. The quantity T = V d 2 V / (d^V) 2 in equation dHJ) 
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TABLE I: Properties of the critical points of the autonomous system (|23[) for the inverse power-law potential V((j>) = Vo4> 
For points M and U the variable v can take any value within the phase space, including v — 0. 

Equilibrium Point x y v Existence 0.$ lu^ q 

M v Always -1 ttn^zl 

dS 10 " 1-1-1 

U ±1 v " 'tto^l 



TABLE II: Eigenvalues of the linearization matrices for the critical points in table [I] 



Equilibrium Point 


X 


y 


V 


Ai A2 


A3 


M 








V 


-3 


37m /2 


dS 





1 





-3 


— 37m 


U 


±1 





V 


6 


37m /2 



is, in general, a function of cf>. The idea behind the 
method in [22J is that T can be written as a function 
of the variable v, and, perhaps, of several constant pa- 
rameters. Indeed, for a wide class of potentials the above 
requirement - T = T(v) -, is fulfilled. Let us introduce a 
new function g(v) — v 2 (T(v) — 1) so that equation (|2Tj) 
can be written in the more compact form: 

v' = -VGxyg(v). (22) 

Equations (fT5|) . (fTl|) . and (f2"2"j) form a three-dimensional 
(closed) autonomous system of ODE: 

x' = (x 2 — l){3a; — Viyv}, 

1 V co R 3(x 2 - j m )y 2 

V = oi 3 7m - V3xyv H }, 

^ VI - a; 

u' = —V3xyg(v), (23) 

that can be safely studied with the help of the stan- 
dard dynamical systems tools [23j]. The function 5(f) 
can be analytically written for a wide variety of poten- 
tials. If g(v) were a polynomial in v (it is the case for 
most quintessential potentials of cosmological interest), 
then, with each root v — v oi of the polynomial equation 
g(v) = (the voi-s include the non-vanishing roots of the 
polynomial as well as the trivial root v = 0), it can be as- 
sociated an equilibrium point of the autonomous system 
of ODE dU). When v = 0, V(<f>) = V , while, if v 0l ^ 0, 
then Vi(4>) oc exp (— t»oi</>). Therefore, equilibrium points 
that arise under the requirement g{v) = 0, are associ- 
ated either with exponential potentials in the DBI field 
4 1 ( v 0i 7^ 0), or with the constant potential (v = 0). 

A rough inspection of the ODE (|2"3"|) reveals that, inde- 
pendent on g(v), for critical pints with y = 0, then, either 



x = ±1, or 1 = 0. Therefore, independent on the func- 
tional form of the potential V(<f>), the kinetic/potential 
energy-scaling solution (x = ±1 => <fi 2 = V(<p)), and 
the matter-dominated solution (x = 0, y = 0, i. e., 
3-ff" = p m ) are always equilibrium points of the au- 
tonomous system of ODE (f2"3")l . Alternatively, for equi- 
librium points with x = (j/ ^ 0), 

2/ = v^yv, y' = -y i y( 1 - 2/ 2 ) > 

so that, necessarily y = ±1, v = 0. This means that for 
polynomials 3(f), for which u = is a root of the polyno- 
mial equation g(v) = 0, the de Sitter solution 3H 2 = Vq 
is an equilibrium point of the autonomous system (|23|) . 

If one started with standard tachyon cosmological 
equations (|11[) instead, then, due to the square root 
y/l — V> 1 one were forced to choose the phase space vari- 
ables [il 

The latter variable is necessary to close the corresponding 
system of ODE (compare with our variable v). Contrary 
to the case with the variable v , it can be shown that it is 
very difficult to write the parameter T = Vd 2 V/ (d v V) 2 
for arbitrary potentials, as a function of the variable 
A. In fact, only power-law potentials, in particular 
V(tp) cx <p~ 2 , can be analytically investigated [l4|. To 
study the dynamics driven by other self-interaction po- 
tentials one has to rely on the obscure concept of "in- 
stantaneous critical points" . 

At this point one recognizes the importance of the 
DBI/tachyon equivalence discussed before. Actually, as 
already shown, the cosmological equations of standard 
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FIG. 2: A plot of the variable v vs <j) for the potential V = 
Vq4>~ x is shown for the chosen values of the free parameters: 
Vo = 1, A = 2 (left-hand panel), while the same plot for the 
potential V = Vo smh~ a (X(j)) is shown in the right-hand panel 
for Vo = 1, A = 1 and a — 3. 



77 = A0-\ (25) 
so that the following asymptotics hold true: 



lim v(d>) 



±oo 



lim v((j>) = 0. 



(26) 



In the left hand panel of Figure[2]a plot of v(<f>) vs is 
shown for the chosen values of free parameters: Vo = 1, 
A = 2. Notice that the range </)£] — oo, 0[ is covered by 
negative values of the variable v, while positive values 
of v > cover the range </> G]0, oo[. In what follows, 
for definiteness, we will restrict ourselves to » > 0, so 
that the tachyon field variable takes values in the interval 
<j> e]0, oo [. In this case the function g(v) in (|22I23|) can 
be written in the following way: 



tachyon cosmology are mathematically equivalent to 
© under the change of variable (fTU|) . Hence, the present 
approach based on the variables defined in (fT2|) . ([20)) (in 
particular the variable v), enables us to further general- 
ize former studies - with the application of the dynamical 
systems tools to tachyon cosmology - to a wider variety 
of self-interaction potentials. 

A drawback of the approach undertaken here can be 
associated with the fact that, for potentials that vanish at 
the minimum - usually correlated with late-time dynam- 
ics - the variable v is undefined so that the corresponding 
behavior can not be properly investigated in the phase 
space spanned by x,y,v. 

The phase space where to look for equilibrium points 
of the system of ODE Ip3[). corresponding to the MTC 
model, can be defined as follows (we take into account 
only expanding universes so that only y > are being 
considered): 

* = {(x, y, v) : -1 < x < 1, < y 4 < 1 - x 2 , v}, (24) 

where it has to be pointed out that the range of the 
variable v depends on the specific kind of self-interaction 
potential considered. Recall that, for constant and ex- 
ponential self-interaction potentials one does not need to 
consider the latter variable, so that, the corresponding 
system of ODE is a two-dimensional one. 

In the next subsections we study particular examples 
were the usefulness of the approach undertaken here is 
illustrated. 



1. The (inverse) power-law potential V((f>) = Vo<f> X 

The inverse power-law potential have been extensively 
studied within standard tachyon field model [3, HH . Ac- 
cording to the definition ([2"U[) of the variable v, for this 
potential one has: 



g(v) = v 2 /X. (27) 

The cosmic dynamics driven by this potential can 
be associated with a 3-dimensional phase space (|24p . 
spanned by the variables x, y, and v, where < v < oo. 

The equilibrium points of the autonomous system of 
ODE ([2"5[) in the phase space ^ defined above, are listed 
in table U while the eigenvalues of the corresponding lin- 
earization matrices are shown in table |TTJ These are non- 
hyperbolic critical points meaning that only limited in- 
formation on their stability properties can be retrieved 
by means of the present linearized analysis. To help us 
understanding the stability properties of these points one 
has to rely on the phase portraits which show the struc- 
ture of the phase space through phase path-probes orig- 
inated by given initial data. 

Existence of the matter-dominated solution (equilib- 
rium point M in table [TJ> , is independent on the value 
of the variable v, meaning that this phase of the cos- 
mic evolution may arise at early-to-intermediate times 
(0 < v < oo), as well as at late time (v — 0). As seen 
from table [Til since in this case the two non vanishing 
eigenvalues of the linearization matrix are of opposite 
sign, the matter-dominated solution is always a saddle 
equilibrium point of (1231) . This solution is inflationary 
whenever 7 m < 2/3. 

Something similar can be said about the equilibrium 
point U in TabUJ This point can be associated also ei- 
ther with early-time, as well as with intermediate-time 
dynamics, due to the fact that v can be any value. How- 
ever, unlike the matter-dominated solution, the critical 
point U could be a past attractor (source critical point, 
usually associated with early-time dynamics) in the phase 
space, since the two non- vanishing eigenvalues of the cor- 
responding linearization matrix are positive. Notwith- 
standing, since the point is a non-hyperbolic one, no 
final statement can be made on its stability properties 
until the corresponding phase portraits are drawn. The 
point U is also a matter-dominated solution, however, 
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TABLE III: Properties of the critical points of the autonomous system (|23|l for the potential V = Vb[sinh(A0)] -a . 

Here v > aX (the constant potential v = corresponds to the particular case where either a = 0, or A = 0), and we introduced 

the following constant parameter: y 2 = (\/36 — q 4 A 4 — a 2 A 2 )/6. 



Eq. Point 


X 


y 


■f 


Existence 




W0 


3 


M 








U 


Always 





-1 


(3 7m - 2)/2 


V 


±1 





D 










(3 7m - 2)/2 


T 


aAy./v^ 


y* 


Q:A 




1 




(3 7m -2-a 2 A 2 y 2 )/2 


MS 


V7m 


y/ 37 m /a A 


aA 


3 7m < ct 2 \ 2 y 2 


37m/a 2 A 2 Vl - 7m 


— 1 + 7m 


(3 7m - 2)/2 



TABLE IV: Eigenvalues of the linearization matrices cor responding to the equilibrium points in Tab lHIl 
Here n = V (48 7m y/l - 7 m/a 2 ^) +4 + 7m (17 7m - 20). 



Eq. Point 


X 


y 




Ai 


A 2 


A3 


M 










-3 





37m /2 


U 


±1 
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3 7 m/2 


T 


aAy»/^ 




a A 


-2a\ 2 y t 


-3 + a 2 A 2 y 2 /2 


-3 7m + a 2 \ 2 y 2 


MS 




V37m/aA 


a A 


—67m /a 


-3[(2 - 7m ) + n]/4 


-3[(2- 7m )-n]/4 



this phase corresponds to the ultra-relativistic case (the 
Lorentz boost 7 — > 00) and, unlike the standard situ- 
ation, it represents scaling between the kinetic and the 
potential energy densities of the tachyon scalar. As be- 
fore, the solution is inflationary whenever 7 m < 2/3. 

The late-time dynamics driven by the potential V = 
Vo<fi~ x is correlated with infinitely large values of the vari- 
able </>—>■ 00 which, in the phase space (|2"4l (0 < v < 00), 
is depicted by the equilibrium point with v = in table 
U (equilibrium point dS). This point corresponds to the 
inflationary de Sitter solution (3H 2 = V). From TabHIl 
it is seen that this equilibrium point could be a late- 
time attractor since the two non- vanishing eigenvalues of 
the corresponding linearization matrix are both negative. 
However, since as already said, this is a non-hyperbolic 
point, only after drawing the corresponding phase por- 
traits one is able to make conclusive statements about its 
stability properties. 

We want to notice that the above results remain true if 
one considered power-law potentials with negative values 
of the constant parameter A. In particular the above 
results can be safely extended to the quadratic potential 

V((j)) (xcf) 2 . 

In FigOJ in order to illustrate the stability proper- 
ties of the asymptotic solutions M, dS, and U, probe 
paths in phase space - trajectories in the phase space 
originated by given initial data - are drawn. As clearly 
seen, these trajectories emerge from the ultra-relativistic 
(matter-dominated) point U and converge towards the 
inflationary de Sitter attractor (point dS) at late times, 




FIG. 3: In the left-hand panel trajectories in phase space for 
different sets initial conditions are drawn for the potential 
V = Vo</!>~ A (Vb = 1, A = 2, 7 m = 1 - background dust), 
while the flux in time of the corresponding system of ODE is 
depicted in the right-hand panel of the figure. Probe-paths in 
the phase space originate at the ultra-relativistic equilibrium 
point U {x = ±1) while at late times approach to the de Sitter 
attractor dS (x = 0, y = 1). 



thus confirming our suspects that U is the past attractor, 
while dS is the future attractor. 



2. The potential V (cf>) = V [swh(\ct>)] a . 

This potential has been formerly studied in Ref.(25j 
as a new cosmological tracker solution for quintessence. 
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According to the definition ([20]) . for this potential one 
gets: 

v = aAcoth(A0), (28) 
from which it follows, in particular, that 

lim v(4>) = oo , lim v((f>) = ±aX. (29) 

0— >0 <fr— >±oo 

In the right-hand panel of Figj2]a plot of i>(</>) vs <f> is 
shown for the chosen values of free parameters: Vo = 1, 
A = 1 and a = 3. Notice that the range of the variable 
d> e] — oo,0[ is covered by — oo < v < — <xX, while the 
range cj> e]0, oo[ is covered by aX < v < oo. In what 
follows, for definiteness we will restrict ourselves to the 
interval aX < v < oo. For the above potential the func- 
tion g(v) defined in (|22I23I) can be written in the following 
way: 




FIG. 4: Phase portrait for the model driven by the potential 
V — Vb[sinh(A0)]~ a (left-hand panel), and the flux in time 
of the system of ODE (|23p for this case (right-hand panel). 
The free parameters chosen are: 7 m = .2, A = 1 and a = 1. 
It is clearly seen that the ultra-relativistic solution U is the 
past attractor while, due to the above choice of parameters, 
the matter-scaling solution MS is the late-time attractor. 



, . v 2 - a 2 X 2 , . 

g(v) = . 30 

a 

The cosmic dynamics driven by V — Vb[sinh(A</>)]~ Q 
can be associated with the 3-dimensional phase space 
(l24l) . where the variable v is constrained to the interval 
aX < v < oo. The equilibrium points of the autonomous 
system of ODE (|2"3"|) in the phase space "J defined above, 
are listed in table IIII1 while the eigenvalues of the corre- 
sponding linearization matrices are shown in Tab II VI 

As for the power-law potential, the existence of the 
matter-dominated solution (equilibrium point M in ta- 
ble IIIII) . is independent of the value of the variable v, 
meaning that this phase of the cosmic evolution may 
arise at early-to-intcrmcdiatc times (aX < v < oo), as 
well as at late times (v = aX). As seen from Tab lIV| 
since in this case the two non vanishing eigenvalues of 
the linearization matrix are of opposite sign, the matter- 
dominated solution is always a saddle equilibrium point 
of ([23)) . The corresponding cosmological solution repre- 
sents decelerating expansion whenever 7 m > 2/3. Un- 
like this, the matter-dominated equilibrium point U in 
Tab lllll can be associated with ultra-relativistic behavior 
(large Lorentz boost). As already said this point repre- 
sents scaling between the potential and the kinetic ener- 
gies of the tachyon field. As it can be seen from the phase 
portraits, it is always the past attractor for any path in 
the phase space of the model. 

Equilibrium points T (the tachyon-dominated solu- 
tion) and MS (the matter-scaling solution) are associ- 
ated with late-time dynamics since, according to (|28[). 
v = a A is correlated with infinitely large values of the 
variable <j). The scalar field-dominated solution T always 
exits and whenever -y m > a 2 A 2 (\/36 + a 4 A 4 - a 2 A 2 )/18 
it is a stable equilibrium point (the late-time attrac- 
tor), otherwise it is a saddle critical point in phase 
space. Whenever the matter-scaling solution MS exists, 



it is a stable equilibrium point (the late-time attractor). 
This solution is always associated with accelerated ex- 
pansion. As in Ref. 14], one has to take caution since 
the critical point MS does not exist if either 7 m = 1 
(matter-dominated era), or 7 m =4/3 (radiation domina- 
tion). This is due to the fact that the existence of the 
matter-scaling solution requires fulfillment of the condi- 
tion < 7,„ < 1. In this sense this solution can not be 
associated with a realistic model of dark energy. 

In the figure |4] the phase portrait for this case is de- 
picted. The above discussed behavior is clearly illus- 
trated by the figure. The free parameters were taken in 
such a way that the MS solution exists (j m = .2, A = 1, 
a = 1). In correspondence, the ultra-relativistic phase 
U is the past attractor, while the matter-scaling solution 
MS is the late-time (inflationary) attractor. 

IV. DBI/TACHYON EQUIVALENCE 

As already shown, there is a mathematical equivalence 
between the DBI model given by the Lagrangian density 
(O, and the tachyon model portrayed by the Lagrangian 
([T]), under the transformation (flQ)) : 



-> V = 




Equivalence under (|10[) implies a residual equivalence be- 
tween magnitudes of physical relevance, in particular: 



ft = Q,<p, q(4>) =q((p) = -(1 + Tp)- (31) 

It is evident that, once the functional form of the self- 
interaction potential V{4>) (or V(<p)) is known, the func- 
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tional relationship tp = ip(4>) (or <fi — 4>{f)) can be ob- 
tained through integration in quadratures, so that one 
is able to transform the potential V(4>) = V(cj)(ip)) —> 
V = V{if) (or V(ip) = V{ip((f>)) -> V = V((f>)). Actu- 
ally, from (fTOj) it follows that 



that the non-negative range of the variable v is associated 
with negative values of the tachyon scalar cf) e] — oo, 0[). 
It remains just to translate the corresponding results so 
that one could discuss their physical implications for the 
tachyon cosmological dynamics. 



(p-tp = 




or, cj, - 0o = / \/vW)dip. (32) 



By using equation ([32]) it can be found that the trans- 
formation (|10[) implies the following transformations be- 
tween modified tachyon potentials V((f>) and standard 
tachyon ones V(ip): 

V{4>) = V e-^ ->■ V(tp) = %(<p - ifo)- 2 , (33) 

where Vq = 4/A 2 , and ipo is an integration constant. For 
the inverse power-law potential one gets that 

V(cb) = V <f>- 2X -> V{fp) = %(v - w)- 2 '\ (34) 

where V = [V /{\ + l) 2 *] 1 /^ 1 , and n = A/(A + 1). 
Additionally, for the sinh-like potential oc sinh _2 (A0), 
one obtains the following equivalence 

V{4>) = l/ o sinir 2 (A0) V(cp) = J (35) 

where V = 1/A 2 and ^ = 1/V A 2 . 

Notice that under the transformation (ITUl) the DBI 
exponential potential is equivalent to the square-law 
tachyon potential oc (p studied in [H, while the ( in- 
verse) power-law DBI potential oc <f>~ 2X is equivalent to a 
(inverse) power-law tachyon potential oc (p~ 2n , not fully 
investigated in the same reference. I. e., the inverse 
power- law potential is not transformed under (|10[) . In 
the later case, the only difference of physical significance 
is in the power of the potential since A— ^n = A/(A+l). 
It is seen that, for positive A > =>• n < 1. 

The above discussed DBI/tachyon equivalence opens 
up the possibility to apply the present approach to inves- 
tigate the dynamics of standard tachyon cosmology for 
self-interaction potentials beyond the square-law poten- 
tial which has been studied in detail in Ref. [l4| . Actually, 
consider, for instance, the exponential self-interaction po- 
tential for the standard tachyon: 

V{ip) = Vbe w . 

By using the relationship (|32l) it can be shown that: 

V(cp) = VoeW -> V{4>) = 

Fortunately this case for the DBI-modified dynamics has 
been already studied in subsection IIIIB II (just replace 
A —> —2 in equations (|23|) and (|27l) . and bear in mind 



V. DISCUSSION 

Thanks to the formal mathematical equivalence un- 
der the transformation (ITU)) , among standard tachyon dy- 
namics depicted by Sen's Lagrangian ((T|) and the dynam- 
ics of a DBI field given by the Lagrangian within action 
([7]). the approach undertaken in this paper enables ap- 
plying the standard tools of the dynamical systems to 
investigate the cosmic dynamics driven by a wide vari- 
ety of (scalar) tachyon self-interaction potentials, with- 
out resorting to such obscure concepts as "instantaneous 
critical points", whose physical relevance is suspicious. 
Actually, if such a mathematical (and dynamical) equiv- 
alence is taken into consideration, the results obtained in 
section III - after applying the linear analysis to study 
the dynamics of the model of ([7]) - can be safely trans- 
lated to the case of the standard tachyon model portrayed 
by the Lagrangian density (TIJ. 

As shown in the former section there is a full equiva- 
lence between inverse power-law DBI potential oc </>~ 2A 
and that of the tachyon oc <p~ 2n , so that, for this kind of 
potential the results displayed in the tables [T] and QTJ for 
the DBI field hold true for the standard tachyon cosmo- 
logical model of [14| , which means in turn, that a detailed 
study of the dynamics driven by this tachyon potential 
is indeed possible. According to the results of section 
III (see Tabs. HJand|TI|, for the potential V(<p) oc <^ 2n , 
whenever < n < 1, one obtains that the de Sitter so- 
lution - point dS in Tab. [T] — is always the late-time 
attractor in the phase space, while the ultra-relativistic 
matter-dominated solution - point U in Tab. [TT] — is the 
past attractor from which the phase paths originate. The 
matter-dominated solution M is always a saddle in the 
phase space. Therefore, the standard tachyon cosmology 
model driven by the inverse power-law potential could be 
a nice scenario to address the late-time cosmic accelera- 
tion. 

The study of the asymptotic properties of the tachyon 
model for the exponential potential is mathematically 
equivalent to the study of the asymptotic properties of 
the DBI cosmological model for the quadratic potential 
V((j>) oc 4> 2 , which is a particular case of the study pre- 
sented in 1111 B II when the constant parameter A is re- 
placed by the particular negative value —2. The only 
think to be changed is the phase space itself if one keeps 
4> S]0, oo [ — in this case, in place of the half of the phase 
space (|2~4"]l corresponding to positive v-s, one has to con- 
sider the complementary half defined by negative v-s in- 
stead — , or one might keep intact the phase space at 
the cost that the tachyon field itself takes values in the 
interval cf> g] — oo, 0[. 
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A remarkable property of the tachyon model studied 
here is that, independent of the particular functional 
form of the self-interaction potential V(4>) considered, 
the matter-dominated solutions M and U - the ultra- 
rclativistic matter-dominated solution, are always equi- 
librium points of the corresponding autonomous system 
of ODE (H3J (see tables HUTU) . A straightforward inspec- 
tion of the equations (I23|) reveals why this happens. Ac- 
tually, a crude inspection of the equations in the system 
of ODE ([25)1 shows that, independent of the functional 
form of the function g(v) and of the value of the vari- 
able v, since for y = the system (|23l) reduces to the 
simplified system of equations: 

x =3x(x 2 -1), y' = Q, v' = 0, 

then, for x — and x = ±1, the corresponding points 
(x, y, v) in phase space: M = (0, 0, v), and U = (±1, 0, v), 
both are equilibrium points of the system of ODE (1231) . 
Since the existence of these points is independent of the 
value of the variable w, both phases of the cosmic evolu- 
tion may arise at early, intermediate, as well as at late 
times. In fact, the point M is always a saddle critical 
point, while U is the past attractor for any path in the 
phase space of the model, otherwise, U is the point in 
phase space from which all of the phase trajectories are 
repelled. 

From the analysis of the equations (|2"3"|) it also arises 
that, in general, for potentials for which the function g{v) 
is a polynomial in v (it happens for most part of known 
quintessential potentials), and the polynomial equation 
g(v) — has non-vanishing roots v = «oi 7^ 0, since 
in this case the system ([23]) reduces to the autonomous 
system of ODE (TT5)) for an exponential potential V oc e A(?i 
(A = vq), the late-time dynamics of the tachyon field 
can be either the scalar field-dominated solution, or the 
matter-scaling phase. This conclusion is quite robust and 
has been formerly stated in [22j in a different context. 

We can use the formal DBI/tachyon mathematical 
equivalence to establish links between previous results 
within the DBI bibliography and the study of equivalent 
tachyon field cases. In reference [2l[, for instance, the au- 
thors study (among other cases) the late-time dynamics 
of a DBI model based on the action @ , with exponential 
potential and brane tension (here we use the notation of 

my- 

V{<j>) = ae-^, }{cj>) = 7 e-"* . 

In the particular case where A + /i = => f • V = 
7 cr, a new class of solutions arise, for which < 7 < 1 
(7 = 7 _1 = Vl — x 2 ), i- e., x 7^ 0, and x ^ ±1. These 
solutions are the matter-scaling solution and the scalar 
field-dominated (kinetic/potential energy-scaling) phase. 
It is interesting noting that, under the field replacement 
(jTUJ) , and / • V — 1, the exponential potential above 
transforms into the (inverse) square tachyon potential 



so that, given that ja = 1, the former results are equiv- 
alent to the results of the dynamical systems study of 
tachyon cosmology with an inverse square potential of 
reference [lij]. Actually, the equilibrium points (c) and 
(dl,d2) in Tab. I of Ref. [l4| are the above mentioned 
matter-scaling, and the scalar (tachyon) field-dominated 
solutions respectively. As long as we know, no such equiv- 
alence between DBI and tachyon dynamics has been re- 
ported before. 

VI. CONCLUSIONS 

In the quest for alternative models of inflation (here 
we include both primordial early-time and late-time in- 
flation), the so called tachyon models and their general- 
ization: the scalar DBI field models, have played an im- 
portant role. Due to complexity of the analysis, in partic- 
ular in connection with difficulties to obtain closed (au- 
tonomous) systems of ordinary differential equations out 
of the corresponding cosmological equations, the study of 
the tachyon dynamics has been performed for a limited 
number of particular tachyon potentials. 

In the present paper we have proposed an approach 
based on an existing mathematical equivalence between 
the dynamics of a DBI field <f> with the following relation- 
ship between the warp factor / of the internal geometry, 
and the DBI potential V: / • V = 1, and the dynamics 
of a standard tachyon field ip driven by the Lagrangian 

under the field redefinition ip — J dej)/ yV(</>). The above 
equivalence allows to translate the results of dynamical 
systems studies of DBI models with the condition / • 
V = 1, for given DBI potentials, to the corresponding 
studies of the equivalent tachyon models. The discussed 
equivalence, in conjunction with the application of an 
approach developed in Ref. 22], which allows to apply 
the standard dynamical systems tools to the study of a 
wide variety of self-interaction potentials, provided that 
the function 

g(v) = v 2 (T-l), 

where 

= J_K y Vd l V 
v ~ V ' {d^Vf ' 

can be written as a polynomial in the new phase space 
variable v, makes possible to investigate tachyon models 
driven by a broad class of potentials. The usefulness 
of our combined approach has been illustrated with the 
detailed study of the following tachyon potentials 

V(v) = V QV - 2n , V(<p) = -^- s , 



where the latter potential is just a particular case of the 
former one for n = 1, ipo = 0. It has been demonstrated 
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that the (p~ 2n potential is equivalent to the inverse power- 
law DBI potential (\r 2X (n = A/(A + 1)), while the in- 
verse square tachyon potential ip~ 2 is equivalent to the 
exponential DBI potential exp— \<fi. The particular in- 
verse square tachyon potential (ip 2 — Vo) _1 i s equivalent 
to the sinh-like potential sinh -2 (\cj>), therefore, differ- 
ent kinds of DBI potentials might be equivalent to the 
same tachyon potential. In spite of the fact that only 
the dynamics driven by the inverse power-law tachyon 
potential ip~ 2n has been studied in details, nonetheless, 
the combined method used here can be applied to other 
kinds of potentials, provided that the function g(v) can 
be written as a polynomial on v for the equivalent DBI 
potential V{<j>). 

Amongst the most interesting results of the present 
study we can list the following: 

• A formal mathematical equivalence between DBI 
models with the following particular relationship: 
/ • V = 1, between the warp factor / and the DBI 
potential V, and tachyon cosmological models orig- 
inated from the Lagrangian (fTJ), under the field re- 
definition ([!()]) ■ has been found. Under the above 
DBI/tachyon mathematical equivalence, results of 
different, otherwise unrelated investigations, are 
fully equivalent, so that, for instance, several re- 
sults of studies of DBI models with exponential 
brane tension and potential [21] , and the results of 
the dynamical systems studies of tachyon cosmo- 
logical models with inverse square potential [l4| . 
are revealed to be equivalent . 

• It is revealed that independent of the func- 
tional form of the potential, the matter-dominated 
solution and the ultra-relativistic (also matter- 
dominated) solution, are always associated with 
equilibrium points in the phase space of the tachyon 
models. The latter is always the past attractor, 
while the former is a saddle critical point. 

• It has been demonstrated that, in general, for DBI 
potentials for which the function g(v) can be writ- 
ten as a polynomial in v (it is the case for most 
quintessential potentials of cosmological interest): 

g(v) = £(« - vo) n , 

n 

the late-time dynamics is associated with either the 
de Sitter solution if at the critical point v = 0, 
is a root of the polynomial equation g(v) = 0, or 
with matter-scaling and/or scalar field-dominated 
solutions if at the critical point, the given root v = 

The present approach can be safely (and moderately 
easy) applied to consider new classes of DBI potentials 
that are equivalent to other classes of tachyon potentials 
beyond the (inverse) power-law one. 
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VII. APPENDIX: DYNAMICAL SYSTEMS 

Here we include brief tips of how to apply the dynam- 
ical systems tools to situations of cosmological interest. 
In order to apply these tools one has to follow the steps 
enumerated below: 

1. To identify the phase space variables that allow 
writing the system of cosmological equations in the 
form of an autonomus system of ordinary differen- 
tial equations (ODE), say: 3 

Xi — (^1 i X2 j • ■ ••En) ■ 

2. With the help of the chosen phase space variables, 
to build an autonomous system of ODE out of the 
original system of cosmological equations (r is the 
time-ordering variable, not necessarily the cosmic 
time): 




— fi( x li x 2, ■■■X n ) ■ 



Notice that the RHS of these equations do not de- 
pend explicitly on r (that is the reason why the 
system is called autonomous). 

3. To identify the phase space spanned by the chosen 
variables (x\,x 2 , ■■■x n ), that is relevant to the cos- 
mological model under study. This amounts, basi- 
cally, to define the range of the phase space vari- 
ables that is appropriate to the problem at hand: 

VP = {(xi, X2, ■■■x n ) : bounds on the Xi-s} . 

4. Finding the equilibrium points of the autonomous 
system of ODE, amounts to solve the following sys- 
tem of algebraic equations on (xi,X2, ■■■x n ): 

fi(xi,x 2 , ...x n ) = . 

5. Next one linearly expands the equations of the au- 
tonomous system of ODE in the neighborhood of 
the equilibrium points pk — Pk(xi, x%, ...x n ), k — 



There can be several different possible choices, however, not all 
of them allow for the minimum possible dimensionality of the 
phase space. 
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l,2,...m: 4 I. e., one replaces X{ — > Xi + e^, where 
are the small (linear) perturbations around the 
equilibrium points. Hence the system of ODE be- 
comes a system of linear equations to determine the 
evolution of the e,;-s: 



where the amplitudes (eo)^ are constants of inte- 
gration. 



3=1 



djj 

dx. 



otherwise, since = fi(p) — 0, then 
|l = ^[M(p)|] e j+ 0(ej), 

3 

where we have introduced the linearization or Ja- 
cobian matrix \M\\ = dfi/dxj. 

6. The next step is to solve the secular equation to de- 
termine the eigenvalues Xi of the linearization ma- 
trix at the given equilibrium point p: 

det\M(p)i-\Uf\=0, 

where [U?] is the unit matrix. 

7. Once the eigenvalues of the linearization around a 
given equilibrium point p have been computed, the 
evolution of the perturbations is given by 

n 

ei(-r) = ^(eo)- exp(Aj-r) , 

4 In general the number of equilibrium points is different from the 
dimension of the phase space: m ^ n. 



If all of the eigenvalues have negative real parts, the per- 
turbations decay with r, i. e., the equilibrium point is 
stable against linear perturabtions. The corresponding 
equilibrium point is said to be a future attractor. If at 
least one of the eigenvalues has positive real part, the 
perturbations grow with r so that these are not stable in 
the direction spanned by the given eigenvalue. Hence the 
point is said to be a saddle. The perturbations around 
a given equilibrium point are unstable, in other words 
the point is a past attractor (a source point in the phase 
space), if all of the eigenvalues have positive real parts. 
Points whose linearization is characterized by complex 
eigenvalues are said to be spiral equilibrium points, and 
are commonly associated with oscillatory behavior of the 
corresponding solution. If at least one of the eigenval- 
ues has a vanishing real part, the equilibrium point is 
said to be non-hyperbolic. In the latter case, in general, 
and unless some of the non-vanishing real parts of the 
eigenvalues are of opposite sign, one can not give conclu- 
sive arguments on the stability of the equilibrium point. 
Other techniques have to be applied. 
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